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ABSTRACT 

We present results of submillimeter observations of a low-mass young 
stellar objects (YSOs) condensation in the Cederblad 110 region of the 
Chamaeleon I dark cloud with Atacama Submillimeter Telescope Experiment. 
Our HCO + (J=4-3) map reveals a dense molecular gas with an extent of ~ 0.1 
pc, which is a complex of two envelopes associated with class I sources CedllO 
IRS4 and IRS11 and a very young object Cha-MMSl. The other two class I 
sources in this region, IRS6 and NIR89, are located outside the clump and have 
no extended HCO + emission. HCO + abundance is calculated to be 2.6 x 10~ 10 
for MMS1 and 3.4 x 10~ 9 for IRS4, which are comparable to the reported value 
for other young sources. Bipolar outflows from IRS4 and IRS6 are detected in 
our 12 CO( J=3-2) map. The outflow from IRS4 seems to collide with Cha-MMSl. 
The outflow has enough momentum to affect gas motion in MMS1, although no 
sign has been detected to indicate that a triggered star formation has occurred. 

Subject headings: stars: formation — stars: low mass — stars: pre-main sequence 
- ISM: clouds — ISM: jets and outflows — stars: individual (CedllO IRS4, Cha- 
MMSl) 
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INTRODUCTION 



1.1. Single Star Formation and Cluster Formation 



A number of past observations and theoretical researches have revealed a rough scenario 
for the formation of isolated, single low-mass stars. That is, a gravitational collapse of 
a molecular cloud core is the startin g point of star formation, then a protostar with a 
thick envelope and bipolar outflows ( iBachilleii Il996l ) is born. The outflows blow off the 
protostell ar envelope, which coul d have an impact on the mass assembling process of the 
star (e.g. lArce fc Goodman! 120021 ). affect the motion of the cloud ar ound the protostar (e.g . 



Takakuwa et al 



200 



^), and cause outflow-triggered star formation ( Yokogawa et al. 2003 ) 



as well as generates turb ulence and limit the efficiency of star formation (INakano et al.lll995l ; 
Matzner fc McKeelboOok 



It is still veiled how to form stars in a group or a cluster (ILada &: Ladall2003l ). although 
most stars are known to be born that way. If stars are forming in a group or cluster, 
impacts of outflows on the environment may be greater than in the formation of isolated 
stars. Therefore, it is important to study interactions between outflows and dense gas in 
cluster-forming regions for comprehensive understandings of star formation processes. 



1.2. Cederblad 110 Region 



The Cederblad (Ced) 110 region (jCederbladl 119461 ) is located in the center of the 



Chamaeleon (Cha) I dark cloud, one o f the low-mass star fo rmation regions in the solar 
vicinity (£>=160pc; [whittet et al1ll997h . Iflaikala et all tood ) identified a C 18 clump in 
the Ced 110 region, whose radius is 0.19 pc and mass is 11.7 M m . This is the most massive 



clump in their samples. Infrared and millimeter continuum observations ( 


Pri 


isti e 


t al. 


1991; 


Reipurth et al. 


1996; 


Zinnecker et al. 


1999: 


Lehtinen et al. 


2001; 


Per si et al. 


2001 


) have re- 



vealed that seven YSOs at different evolutionary stages and one very young source are 
clustered within ~0.2 pc. In Table [T] we list the objects with their position, IR class, T ho \ 
f lMvers fc Laddlll993h . and L hoh 



The young est object in th i s reg ion, Cha-MMSl, was detected by 1.3 m m continuum 
observations by lReipurth et al.l (119961 ). At this position. iLehtinen et al.l (120011 ) identified the 
far-IR object Ced 110 IRS10 by observations with the Infrared Space Observatory (ISO), 
while no corresponding source has b een detected with nea r-IR observ ations. Thus, MMS1 is 
considered to be a "real protostar" (ILehtinen et al.ll200ll ). However, ILehtinen et al.l (120031 ) 
could not detect 3 or 6 cm radio continuum emission. This result indicates MMS1 does not 
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have any jets or outflows. They noted that MMS1 is still at a prestellar stage. 



CedllO IRS4, a class I object, is associated with a reflection neb ula that is extended in 



a north-south di re ction and is div i ded b y a thick edge-on dust disk (IZinnecker et al.lll999 



Persi et al.ll200ll ). iHenning et al.l (119931 ) detected 1.3 mm radio continuum emission from 
this object, while extended emission was not detected (IReipurth et al.l 1 1 9961 ) . 



Molecular outflow was detected in this re gion by 12 CO (J= 1- ) observations wit h the 



Swedish-ESO Submillimeter Telescope (SEST; iMattila et al. 



19891 : IPrusti et al.lll99lf ). but 



the identification of its driving source has been controversia l. Originally, IRS4 was thought 
to be driving this outflow. However, IReipurth et al.l (119961 ) detected MMS1 and suggested 
that MMS1 is the drivin g source of the outflo w as we ll as the nearby Herb ig-Haro (HH) 
objects HH49/50 (see also iKontinen et al.ll20000. Wh ile iLehtinen et al.l (120031 ) revealed that 



MMS1 is associated with no jets, iBally et al.l (120061 ) suggested that a parsec-scale outflow 



is ejected from MMS1 based on the positional relationship among the HH objects and the 
YSOs in their large area shock survey. Lack of high-resolution mapping observations of the 
outflow (e.g. the mapping grid was 1' in Mattila et al. 1989) have made it difficult to resolve 
this issue. 



Cedll O IRS11 was first det ected by IZinnecker et al.l (119991 ) as IRS 4E in their near- 
IR image. ILehtinen et al.l (120011 ) tentatively detec ted a far-IR source from ISO data, and 
identified LRS11 as ISO-Chal 86 JPersi et aljboooh. The coordinates for LRS11 in Tabled] 
are based on the position for Chal 86 in lPersi et al.l (120011 ). 



In this paper we focus on the following three points; (1) the distribution of dense gas 
in the YSO clustering region, (2) the influence of outflows from older protostars on newer 
ones and the surrounding gas, and (3) th e " blue asymmetry" line profile tow ard YSOs which 
indicates the existence of infalling gas (ILeung fc Brownl Il977l ; IZhoul |l992j) . To investigate 
them we observed CO(J=3-2) and HCO + (J=4-3) lines. The observations of the CO line 
can reveal the distribution of molecular outflows with high spatial resolution and enable us 
to distinguish the driving source of out flows. The H CO + line has a high critical density of 



9.7xl0 6 cm 3 at the J=4-3 transition fevans 1999 ). and its observations can probe dense 



gas and also can detect gas infall in protostellar systems through detection of the blue 
asymmetry in line profiles. We also made observations of the optically thin H 13 CO + line 
toward CedllO IRS4 and Cha-MMSl to obtain precise physical parameters. In addition, 
we observed CH 3 OH.(J=7k-6k), SiO(J=8-7), and SO (J=8 g -7 8 ) lines, whose upper energy 
levels^] are 70.5, 75.0, and 87.5 K, respectively. These high energy levels made these lines 
useful diagnostic tools of shocked regions. 



1 The data are taken from the Leiden Atomic and Molecular Database (LAMDA). 
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We introduce the detail of our observations in §2, and the results of the observations 
are described in §3. In §4 we discuss molecular abundance, gas infall toward protostars, 
kinematics of outflows, and the interaction between MMS1 and the outflow from IRS4. Our 
conclusions are summarized in 55. 



2. OBSERVATIONS 

The CO( J— 3-2) and HCO + ( J =4-3) observations were perform ed in 2004 Novembe r with 



the 10 m Atacama Submillimeter Telescope Experiment (ASTE, lEzawa et al.ll2004j ). Ob- 
servations were remotely made from the ASTE operation rooms at San Pedro de Atacama, 
Chile and National Astronomical Observatory of Japan (NAOJ) Mitaka campus, Japan, 



using the network observation system N-COSMOS3 developed by NAOJ (IKamazaki et al. 



20051 ). We used a 345 GHz (0.8 mm) band double side band (DSB) superconductor-insulator- 
superconductor mixer receiver using position switching. Taking advantage of the DSB obser- 
vation, we simultaneously received 12 CO(J=3-2; 345.796 GHz) in the lower side band (LSB) 
and HCO + (J=4-3; 356.734 GHz) in the upper side band (USB). With the same system, 
we observed H 13 CO+( J=4-3; 346.999 GHz) line in 2005 June and the CH 3 OH( J=7^-6^; 
338.344 GHz), SO (J=8 9 -7 8 ; 346.529 GHz), and SiO(J=8-7; 347.331 GHz) lines in 2006 
December. 

We used the XF digital spectrometer with a bandwidth and spectral resolution set to 
128 MHz and 125 kHz, respectively, except for CH3OH line observations made at 512 MHz 
bandwidth and 500 kHz resolution. The velocity resolution is 0.11 km s _1 and 0.43 km s _1 
in 345 GHz. The half-power beamwidth (HPBW) was 22" at 345 GHz. 

In CO and HCO + observa tions, the mapping g rid spacing was 20" (0.016 pc at the 



distance of Cha I). According to lMizuno et al. 



(119941 ) . the typical size of a star-forming core 



is about 0.03 pc, thus the angular resolution of our observations is high enough to resolve 
this size of a core. With this grid spacing, we obtained spectra toward 136 points which 
resulted in 4' x 3' maps centered at IRS4. The observations were made by position switching 
every 20 or 30 s between the source and off positions. The off positions were (a, <5)(J2000) 
= (ll h 24 m 2r!2, -76°57'39"1) and (ll h 42 m 9:7, -78°03'44"4). These positions were checked 
to be free of emission. 

We observed H 13 CO+ lines toward CedllO IRS4 and Cha-MMSl. The CH 3 OH line 
was observed along the line from IRS4 to MMS1 with 20" grid spacing. SO an d SiO lines 



were observed simultaneously in opposite sidebands toward the iKontinen et al.l (120001 ) SO 
northern peak, (a, 5)(J2000) = (ll h 06 m 34!7, -77°22 / 50"7), and the boundary point between 
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IRS4 outflow and MMS1, (a, 5) (J2000) = (ll h 06 m 34!7, -77°23'24"8). The observed point 
for the shock-tracer molecules are indicated on Fig{TJ 

DSB system noise temperature at an elevation of 35° varied from T sys = 300 to 600 
K. For CO and HCO + lines, on-source integration time for each point ranged from 1 to 5 
minutes, and the typical rms noise was 0.3 and 0.15 K, respectively. We made additional 
integration toward the YSOs for improving signal-to-noise ratio in the HCO + line, reaching 
a rms noise of 0.1 K. For other lines, we made the rms noise level go down to 0.05 K for 
H 13 CO + , 0.1 K for CH 3 OH, 0.15 K for SiO, and 0.2 K for SO. Despite a careful inspection, 
no emissions from shock-tracing molecules such as CH3OH, SiO, and SO were detected in 
our observations. All temperatures are reported in terms of the Tmb scale. The main beam 
efficiency was i] MB ~ 0.6. 

The telescope pointing was regularly checked by observing Jupiter, Saturn, or IRC+10216 
every 1 or 2 hr during observations, and the error was estimated to be a bout 5" peak to 



peak . Intensity calibration was carried out by the chopper-wheel method (IKutner fc Ulich 



19811 ) . The absol ute intensity is ca librated by assuming that Tmb of IRC+10216 for CO 



(J=3-2) is 32.5 K (IWang et al.lll994l ). The daily variability was checked by observing L1551 
IRS5 for HCO+ and by observing Ori KL for H 13 CO+, CH 3 OH, SO, and SiO. The sideband 
ratio was assumed t o be unity. The zenith 220 GHz opacity varied from r ~ 0.05 to ~ 0.1. 



Matsuo et al.l (119981 ) reported the atmospheric transmission is 0.86 and 0.83 in 345 and 356 
GHz, respectively, when the zenith 220 GHz opacity was 0.04. In this paper, we assume the 
atmospheric transmissions in the two side-bands are identical. 

We reduced the data using the NEWSTAR software developed in the Nobeyama Radio 
Observatory. The baseline subtraction was done by a straight line in most cases. 



3. RESULTS 
3.1. Overview 

Figures [T] and [2] show the observed maps of the CO(J=3-2) emission. The CO(J=3-2) 
emission was detected at all points within the 4' x 3' map with a varying degree of intensity. 
Around IRS4 and IRS6, we detected CO high- velocity components most probably caused by 
outflows. The details of the outflows are described in §3.3 and §4.3. 

Figure [3] shows the integrated intensity map of the HCO + (J=4-3) line in the velocity 
range of 2.0 km s _1 < Vlsr < 7.0 km s -1 . A HCO + clump of dense gas is seen in FigJHJ This 
clump is an aggregate of three independent envelopes associated with CedllO IRS4, IRS11, 
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and Cha-MMSl. The envelope associated with IRS 11 seems to be apart from the component 
associated with IRS4 in the Vlsr = 4.78 and 5.00 km s _1 panels of FigJH MMS1 is also 
distinguishable in the Vlsr = 3.92 and 4.13 km s -1 panels. Figure [5] and [H] are the position- 
velocity diagrams of CO and HCO + lines along line A and B in FigJHl respectively. In FigJ5] 
the components associated with IRS4 and IRS 11 are marginally resolved with the viewpoint 
of line width, that is, the emission associated with IRS4 exhibits a distinctly larger linewidth 
than that for IRS11. The HCO+ clump covers 140" x 110" (0.11 pc x 0.085 pc) above the 
3 a noise level. The extension of the clump with a position angle ~ 45 deg corresponds well 
to the 200 /im ridge emission observed with ISO (jLehtinen et al.ll200ll ). 



3.2. CO and HCO+ Line Profiles 

The CO(J=3-2), HCO+(J=4-3), and H 13 CO+( J=4-3) line profiles obtained at YSO 
positions are shown in FigJTl The CO line profile from ambient gas is sampled at the 
boundary of the map and shown in FigJ7]a by a dotted line. 

Almost all CO line profiles have a dip at Vlsr ~ 4.3 km s _1 . The dip also appears in 
some HCO + line profiles at the same velocity. The velo city of the dip is iden tical to that of 



C 18 clump number 3 identified in the CedllO region bv lHaikala et al.l (|2005h. and al s o give s 



close agreement with the velocities of optically thin lines observed by lKontinen et al.l (120001 ) . 



These correspondences mean that this dip is caused by self-absorption by less excited gas, 
not the result of two emission components with different velocities along the same line of 
sight. Since the observed region corresponds to one C 18 clump, the outer gas of the clump 
in the lower excitation level would make the CO absorption. The dip in the HCO + (J=4-3) 
line indicates a significant population of the J = 3 level that requires a moderately high 
density, but with a lower J=4-3 excitation temperature than in the center of the clump. 
This is a reflection of the density gradient in the clump. We can find these dips in Fig. [5] 
and [6j This dip prevents the channel map (Figf4j) from representing the actual distribution 
of interstellar matter in the dip velocity. The indentations of contours located at the south 
of IRS4 and a lack of a peak around MMS1 in the panels of Vlsr — 4.35 and 4.57 km s _1 of 
FigHJmay be due to the self- absorption. It is conceivable that other reasons such as depletion 
and chemistry in the dense gas cause this lack of a peak. The issue of molecular depletion 
is discussed in §4.1. 

The parameters of the lines detected at YSO positions are summarized in Tables [2], [3] 
and HI We fitted the spectral data of HCO + and H 13 CO + with a Gaussian profile. The 
error estimates of the values in the tables come from the fit. In deriving the line parameters, 
the velocity ranges of outflows and self-absorption features in the lines are eliminated in the 
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fitting processes. Although this may have some defectiveness, this is the optimal handling 
of our data. 



3.3. Individual Sources 



3.3.1. CedllO IRS4 and the Outflow 

Bipolar outflow from IRS4 is detected in our CO(J=3-2) observ ations. This outflow 
had been also detected by lMattila et all Jl989h and lPrusti et al.l Jl99lh with their CO(J=l- 
0) observations. We compared the line profiles toward IRS4 with that of the ambient gas 
obtained at the edge of the map in (a, 5)(J2000) = (ll h 06 m 24!3, -77° 21' 1&!7). Judging from 
this comparison shown in FigJ7]a, the line profile in the velocity ranges of Vlsr < 3.0 km 
s _1 and Vlsr >6.3 km s -1 are not affected by the ambient gas, and thus we identified these 
ranges as outflow components. 

The distribution of the outflow components from IRS4 are shown in FigJT] (integrated 
intensity map) and Figj2] (channel map). The blueshifted and redshifted outflows show a 
different aspect. The spatial extent of the blueshifted outflow is considerably larger than that 
of the redshifted outflow. Furthermore, the highest velocity component of the blueshifted 
outflow lies midway between IRS4 and MMS1, wh ile the velocity reaches its peak around 
IRS4 in the redshifted flow. The near-IR images (IZinnecker et al.l Il999l : iPersi et al.l 1200 ll ) 
indicate that IRS4 is seen nearly edge-on. The observed distribution of the blueshifted com- 
ponent of the outflow extending to the south suggests that the outflow axis slightly inclines 
with the southern side closer than the northern side. This configuration is in agreement with 
the previous remarks on the inclination of IRS4 based on the color gradient of the nebulae 
(IZinnecker et al.lll999l ). They found that the northern part of the reflection nebula is redder 
than the southern part, which suggests that the northern part is tilted away from us. Our 
integrated intensity map (FigHJ) d oes not show the r e dshifted componen t in the south of 



MMS1 reported in the literature (IMattila et al.l Il989l ; iPrusti et al. 



199l|) 



In the channel 

map (Figj2]), however, one can see the extended redshifted emission at Vlsr = 6.52 km s _1 , 
and this could be its counterpart. 

Observations of the H 13 CO + ( J=4-3) line were also performed toward IRS4 (FigJT^) with 
a 43 niK noise level. The FWHM of this line, AV = 1.72 ± 0.29 km s -1 , corresponds to that 
of the HCO + line within the margin of error. This supports the detection of the line. With 
this result, we discuss molecular abundance in 54.1. 
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3.3.2. Cha-MMSl 



In Figd?) we show the CO, HCO+, and H 13 CO+ line profiles toward Cha-MMSl. The 
dip in the HCO + line is deeper than that of IRS4. This dee per dip shows that MMS1 is 
deeply embedded in the densest part of C 18 clump number 3 (IHaikala et al.ll2005l ). because 
MMS1 lies nearer to the position of the peak intensity of C 18 than IRS4. 

We calculated the HCO + "original" integrated intensity with a Gaussian fit using the 
line profiles in the velocity ranges which are free from self-absorption. These reconstructed 
" absorption- free" contours are shown with dashed linea in FigOJ The summit ridge of the 
contours extends south and reaches the position of MMS1. This result again supports the 
idea that MMS1 is deeply embedded and its outer envelope makes this absorption feature. 

Around MMS1, there are no conspicuous wide lines (Figj6j). This shows that no outflow 
is associated with MMS1, which is cons istent with the lack of centimeter continuum emission 
from ionized jets (ILehtinen et al.l 120031 ) . Thus, MMS1 is rejected as a driving source of the 
outflows seen in CO (J = 1 — 0) maps (IMattila et al.l 11989c iPrusti et al.l 1 19911 ) as well as a 
launcher of HH49/50. 



There are some literatu r es which summarized properties of clas s objects (IGregersen et al 



19971 ; iMardones et al.l 119971 ; IGregersen et al.l |2000| ; iFroebrichl 120051 ) . We performed a litera- 
ture survey about three aspects of the class objects: centimeter continuum emission, CO 
outflows, and HH objects. All class sources that appeared in the four papers had at least 
one aspect out of the three. On the other hand, MMS1 has no outflow, centimeter continuum 
emission, or HH objects. Our results show no evidence of protostar activity in MMS1. 

The H 13 CO+(J=4-3) line was detected toward MMS1. The width of the line, 0.36 km 
s _1 , is smaller than that of the same line toward IRS4. This optically thin line could reflect 
the cold innermost region of MMS1; thus, the narrow line shows there are no protostars 
formed yet in MMS1. 

The virial mass of MMS1 is derived from the H 13 CO + line width. Assuming the diameter 
of MMS1 is 40" = 0.03 pc jReipurth et alill996h . the viri al mass is calcu l ated t o be 0.49 



Mq. This is comparable to the mass (0.45 M Q ) derived by ILehtinen et al.l (120011 ) from 1.3 
mm continuum observation. 



3.3.3. CedllO IRS 6 



IRS6 is a binary source (IRS6a and I RS6b), and the near-IR luminosity of IRS6a is 
a factor of 7 higher than the companion (jPersi et al.l 120011 ) . This fact indicates that the 
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features obtained toward the binary system mainly derive from IRS6a. 



Toward IRS6, we found a weak outflow component in the CO line (Fig{TJ [SJ and 0c) 
which has not been detected in past radio observations. The redshifted component extends 
north, while the blueshifted component is tightly associated with IRS6. The kinematics of 
the outflow is described in §4.2. 

We marginally detected a weak HCO + line toward IRS6 (Fig{7]c). In FigJH and 0, the 
IRS6 system does not seem to be surround e d by an extended HCO + envelope. This result 
supports previous study by iHenning et al.l (119931 ) that IRS6 has no massive circumstellar 
envelope. The small envelope and high Tboi of 270 K suggest IRS6 is a relatively evolved 
protostar. 



3.3.4. CedllO IRS 11 

Toward IRS11, no high-velocity components were detected by either CO(J=3-2) or 
HCO + (J=4-3) lines. A small dense gas envelope associated with this object is detected with 
the HCO + line (FigOHU and [5]). The radius of the envelope appearing in Figj5]is about 20", 
which corresponds to 1.6 xlO -2 pc at the distance of the Cha I dark cloud. 



3.3.5. NIR89 and Class II /III Sources 



Though NIR89 is classified as a potential class I (IPersi et al.l l200ll ). no CO(J=3-2) 
outflow was detected (Figf7|e). In addition, no HCO + (J=4-3) emission was detected toward 
the source. Because this obj ect is thought of as a "good candidate young brown dwarf with 
a stellar mass of 0.02 M & (IPersi et al.ll200ll ). the amount of the envelope associated with 
NIR89 should be too small to be detected by our HCO + observations. 

No CO outflow components or HCO + emission were detected toward the class II /III 
objects, ISO-Cha I 97, ISO-ChallOl, and CedllO IRS2 fFigfUf - h). This lack of outflow is 
common among T Tauri stars in general. For HCO + emission, Ivan Zadelhoff et al.l (1200 ll ) 
detected the line toward classical T Tauri star LkCa 15 in the Taurus region with a peak T m ^ 
of 0.14 K using the James Clerk Maxwell Telescope. Although the line intensity could vary 
depending on the inclination of the circumstellar disks, in the case of class II /III objects, 
it is still plausible that we could not detect the HCO + line toward those samples in our 
observations with 0.15 K sensitivity (lcr). 
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DISCUSSION 



4.1. Molecular Abundance in Envelopes 

In dense, cold gas envel opes of protostars and prestellar cores, molecules often freeze 
out onto dust grains (see e.g. IjQrgensen et al.1 120051 ) . On the ot her hand, the abunda nce of 
some species like HCO + are known to be enhanced in outflows (IRawlings et al.l I2004T ) . 



We examined the abundance of HCO + toward MMS1 and IR S4 by means of our 
H 13 CO + ( J = 4 — 3) observations and 1.3 mm continuum observations (jHenning et al.lll993l ; 



Reipurth et al.lll996l ). Assuming that the dust emission is optically thin, the column density 
of the H2 molecule, N(H.2), is calculated with 



iV(H 2 



S 



beam 



X 



where S^ eam is the 1.3 mm flux density per beam, f2beam is the beam solid angle, /1 = 2.3 is 
the mean molecular weight, mn is the mass of the H atom, Ki 3mm is the mass absorption 
coefficient per gram, and B U (T D ) is the Planck function at the dust temperature T D . We 
employe d ^.^ mm = 0-01 cm 2 g -1 for IRS4, which is the recommended value for very dense 



regions (jOssenkopf &: Hennind Il994l ). Meanwhile, we applied the value of 0.005 cm 2 g _1 
for MMS1, which is the adopted value for prestellar core L1689B (jAndre et al.lll996l ). The 
substi tuted values of >S',^ eam are 370 mJy for MMS1 (IReipurth et al.lll996l ) and 101 mJy for 
IRS4 (jHenning et al.l 119931 ). We employed a dust tem perature of 20 and 21 K for MMS1 
and IRS4, respectively, based on lLehtinen et al.1 (120011 ). The resulting H 2 column density is 
2.2 x 10 23 cm" 2 for MMS1 and 2.8 x 10 22 cm" 2 for IRS4. 

We also calculated the column density of H 13 CO + under LTE conditions and the as- 
sumption that the excitation temperature equals the dust temperature. Comparing iV(H 2 ) 
obtained from dust observations with iV(H 13 CO + ), we estimated the abundance of HCO + 
to be X(HCO+) = 2.6 x 10" 10 for MMS1 and 3.4 x 10~ 9 for IRS4. We set the isotropy ratio 
of [H 13 CO + ]/[HCO + ] = 70. The results are summarized in Table 



Jprgensen et al.l (120041 ) reported average X(HC O + ) = 8.0 x 10 10 for prestellar cores and 



1.1 x 10~ 8 for class I sources. IjQrgensen et al.l (120051 ) pointed out that the molecular depletion 
in envelopes of class I objects is les s significant because of heating by their central sources. 
Furthermore, IRawlings et al.1 (120041 ) reported that the abundance of HCO + is enhanced in 
molecular outflows. Our results are not in contradiction with abundances reported in similar 
objects. However, it is not appropriate to investigate the results in more detail, because the 
mass absorption coefficient k and the LTE approximation have large uncertainties, k ha s 
different values depending on dust size, shape, structure, and so on (jHenning et al.lll995l ). 
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In addition, iMennella et al.l (119981 ) show k has a dependence on temperature. 



4.2. Gas Infall toward the Protostars 



We show HCO + (J=4-3) line profiles toward four class I sources in Fig. [7]a, c, d and e. 
None of the profile s show the blue asymmetry which is one of the indicators of infalling gas 
toward protostars (ILeung & Brown! 119771 ; IZhoulll992l ). In a simplified model of an infalling 
envelope, the blueshifted line component is stronger than the redshifted component. This is 
because the redshifted component is absorbed by the cooler infalling material in the near half 
of the envelope along the line of sight. Do our results mean these objects have no infalling 
gas? 

As mentioned in the previous sections, IRS4 and IRS6 obviousl y have outflows. This is 
clear evidence of the existence of infalling gas toward the protostars. iGregersen et al.l (119971 . 
2000|) have revealed that blue asymmetry line profiles do not appear in about half of class 
and I sources, although they are thought to be in the active accretion phase. 

Some reasons are suggested why blue asymmetry does not appear in all sources. One 
great problem is confusion with outflows. If a redshifted outflow exists in the telescope beam, 
the redshifted line components absorbed by the infalling envelope would be compensated. It 
is also difficult to separate the influence of other velocity components such as circumstellar 
disks and ambient gas in observations with a large beam of single dish radio telescopes. In 
our observations, as an example, a wing component is detected in the CO line profile toward 
IRS4. This profile obviously shows that our 22" beam covers the outflowing component. In 
general, there are some possibilities to confuse not only outflow, but core rotation and other 
large-scale physical structures with infall, all of which affect obtained line profiles. Thus, 
it is difficult to investigate infall motion in targeted observations toward protostars without 
high-resolution spatial information about the envelope. 

The other problem is optical depth. Choosing molecular lines with appropriate optical 
depth is necessary to detect infall motion of the inner portion of the envelope in general, not 
only for CedllO. When a protostar is highly embedded in a dense clump and an observed 
line becomes optically too thick, the line cannot convey information of the inner infalling 
region any more and reflects only properties of foreground components. 
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4.3. Kinematics of Outflows 

We c alculated th e phys ical parameters of the outflows from IRS4 according to the pro- 
cedures of IChoi et al.l (119931 ) . In these procedures, the CO column density for each channel 
in units of cm -2 , Np°, and the mass per channel, Mj, are calculated from 



N co 



1.10 x 10 



15 T R3 _ 2 A^ 



T~32 



D(n,T K ) 1 - exp(-r 32 ) 



Mi = /im H2 <i 2 fiiV; 



CO 



[H] [C] 
[C] [CO] ' 



where 



(2) 
(3) 



D{n,T K ) = f 2 [J v {T ex ) - J u (T hk )][l - exp(-16.597/T ex )], 



(4) 



Tr3-2 is the line intensity of the CO(J=3-2) line, AV^ is the channel width in km s _1 , / 2 is 
the fraction of CO molecules in the J = 2 state, d = 160 pc is the distance to the Cha I cloud, 
Q is the solid angle subtended by emission, [H]/[C]= 2.5 x 10 3 , [C]/[CO] = 8, and r 32 is the 
optical depth of the li ne. We calculated t 3 ? ~ 0.4 comparing with CO(J=3-2) and CO(J=l- 
0) line-wing intensity (IPrusti et al.lll99ll ). In this calculation, we assumed that the excitation 
temperatures T BY of the two line s are equal to 20 K. According to large velocity gradiant 
simulations by IChoi et al.l (119931 ). D(n, Tk) does not vary so much within the reasonable 
density and temp erature range. Therefore, we set the same value of -D(n,Tk) = 1.5 as in 
Choi et al.l (I1993I ). The mass M, momentum P, and kinetic energy E k can be estimated 
from 



m = >:,,u„ 

P = EMM- V \, 
E K = X-MiiV-Vo) 



(5) 
(6) 

(7) 



where V is the system velocity. We also calculated the dynamical timescale t d = R/V c h 
where the characteristic velocity V c h = P/M, mass loss rate M = M/t&, and average driving 
force F = P/t d . 

The results are sho wn in Table El We adopted an incl ination of the outflow axis from the 
line of sight of % = 72° (IPontoppidan &: Dullemondll2005l ). They derived this value through 
comparing near-IR images with three-dimensional Monte Carlo radiative transfer code of a 
disk-shadow proj ection model. This va l ue is consistent w ith the fact that the high-resolution 



infrared images (IZinnecker et al.l Il999l ; iPersi et al.l l200ll ) clearly show the edge-on disklike 
silhouette of IRS4. 
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W e compared the derived values with those in a statistical study of outflows (IWu et al. 
2004 ) in which they show a correlation between the mass, Lboi and F for 391 samples. The 



mass of IRS4's blueshifted o utflow is smaller by about 1 order of magnitude than the average 
of the same Lboi sources in IWu et al.l (120041 ). but it is still within the range of variation of 



the data. 

In contrast, the derived mass loss rate is rather large in IRS4 despite its small size. 
Generally, the mass loss rate by jet s is about 1% and that by outflows islO % of the mass 
accretion rate (IHartigan et al.lll995l ). The derived value of 1.1 xlO _6 M yr -1 is comparable 
to the typical accretion rate for protostars, 10~ 6 M(p yr" 1 , d erived by the similarity solution 
for thermally supported clouds with 10 K (jHartmann!ll998l ). 



HH49/50 is a group of knot s located 10.5 / southwest o 
dial velocity and proper motion ( Schwartz fc Dopita 198ol; 



knots move toward the southwest and away from us. iMattila et al.l ( 119891 ) and iPrusti et al. 



the CedllO region. Studies of ra- 



Schwartz et al 



19841 ) revealed the 



( 119911 ) identified CedllO IRS4 as the driving source of these HH objects judging from their 
CO outflow map, because they found a diffuse redshifted component that extends in the 
southwest of IRS4 and the alignment of knots is almost parallel to the direction to IRS4. 
In our observations, we detected two components that extend to the south of IRS4. One is 
the prominent blueshifted outflow which has an opposite line-of-sight direction of movement 
from HH objects. The other is the diffuse redshifted component that appeared in FigfSJ 
This complicated morphology makes it difficult to discuss the relation between IRS4 and 
HH49/50. To resolve this issue, one should enlarge CO observations toward the southwest 
to know the larger structure of the outflow. 

In Figj5]and[6j the HCO + line toward IRS4 indicates a wide profi le as well as a CO line . 
Infalling and outflowing gas seem to be the causes of this wide line (IGregersen et al.lll997l ; 
Masunaga fc Inutsukall2000l ). and perhaps, both the effects are blended. In both position- 
velocity diagrams, the blueshifted component of the HCO + line is broader than redshifted 
component. This is identical to the tendency of CO outflow, which indicates the wide HCO + 
line is mainly caused by the outflow. As in the HCO + channel map (FigJ3|), the blueshifted 
component extends to the west of IRS4, whereas redshifted component is elongated to the 
north. This aspect again corresponds with the morphology of the CO outflow (FigfSJ). 

Another method to use to determine what makes the wide line profile is to compare the 
velocity of infalling and outflowing gas. We calculated the free fall velocity of i nf ailing gas, 
fff = .31 - 0.56 km s _1 , assuming the mass range of 0.10 - 0.32 M & derived by lPersi et al. 



( 120011 ). and a diameter of the envelope of 1.7 x 10 2 pc. These are smaller than the line 



width, which indicates the outflow has a larger contribution to the wide line profiles. 
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Toward IRS6, we derived some physical parameters of the outflow in the same procedure 
as for IRS4 and show the result in Table. [6j In the case of IRS6, we could not correct the 
effect of the inclination of the outflow axis because of a lack of information. Besides, there 
have been no CO(J=1-0) observations toward IRS6 in previous years, so we could not derive 
T32. Here we use the same T32 as toward IRS4. Although there are ambiguities of inclination 
and optical depth, the outflow from IRS6 is more than 1 order of magnitude smaller than 
that from IRS4. This is consistent with our estimation of its evolved status with a small 
envelope. This interpretation is not consistent with the young age of the outflow; however, 
when the outflow is nearly in pole-on configuration, the dynamical timescale of the outflow 
will be calculated to be smaller than the true value. We need to know the inclination angle 
of this source to correct this mismatch. 



4.4. Interaction between Cha-MMSl and the Outflow from IRS4 

Figure [T] shows the blueshifted component of the outflow from the IRS4 has negative 
spatial correlation with the continuum emission from Cha-MMSl. The outflow seems to 
change its direction at the edge of the MMS1. In Figj2j the outflow component with the 
highest velocity is distributed at the midpoint between IRS4 and MMS1. In addition, the 
velocity of the outflow suddenly declines in front of MMS1 (FigJBJ), i.e., the flow seems to 
be decelerated by MMS1. Do these characteristics indicate that MMS1 interacts with the 
outflow from IRS4? 

There were no emissions detected from shock tracer molecules in submillimeter wave- 
length in our observ ations. However, considering IRS4's low luminosity (Lboi = 



Lehtinen et al.ll200ll ). weak outflow properties, and high upper energy levels of > 70 K, it is 
conceivable that the outflow-core interaction could not increase the temperature enough 
to excite these molecules. We could not rule out that a weak interaction does occur. 



Kontinen et al.l (120001 ) pointed out that outflows and radiation from the neighboring young 
stars have prospects for influencing the chemical processes in MMS1. Here we considered 
the impact of the interaction from following two view-points. 

First we investigate momentum. Here we compare the momentum brought by the 
outflow with gas momentum in MMS1. When an outflow with density p and velocity v 
is colliding with a surface of radius r for a time duration t, we can obtain the mass M co \ 
colliding in a unit of time and the momentum P ou tflow as, 

-Poutflow = Mcoivt = nr 2 v 2 pt. (8) 



Now we estimate the density of the outflow. For simplicity, we assume the outflow has a cone 
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shape with a base diameter of 40" (= 3.1 xlO -2 pc) and a height of 60" (= 4.7 xl0~ 2 pc) 
judging from FigJTJ Calculating the density p of the outflow from this assumption with the 
derived mass, we obtain p = 6.5 x 10~ 20 g cm" 3 . Assigning this density, radius r = 1.6 x 10~ 2 
pc and t which equals the dynamical timescale of the flow (1.0 x 10 4 yr) to equation (JHJ), 
we derived a momentum of 0.88 M Q km s _1 . Meanwhile, the momentum of gas in MMS1 
is calculated to be Pmmsi — 0.15 M & km s" 1 from the mass of 0.45 M Q and the velocity 
of 0.34 km s" 1 derived from the H 13 CO + line profile, which is a combination of turbulence 
and thermal width, removing the contribution of the spectral resolution. As a result, P ou tflow 
is much larger than Pmmsi ; which means the motion of the gas in MMS1 is easily affected 
when the outflow collide with the core. 



Although the MMS1 contours in lReipurth et al.l (119961 ) have an asymmetry, it is difficult 
to confirm the collisional impact on the geometry of MM SI with their map. But in our HCO + 
position- velocity diagram (Figj6j), we could find the sign of the kinematic motion of the outer 
gas of MMS1 caused by the collision. The dip in the northern side of MMS1 is slightly (~ 
0.3 km s _1 ) blueshifted, which results in the velocity gradient on the map. This gradient is 
consistent with the IRS4 outflow direction. Precise calculations of the momentum for this 
gradient is difficult with our data set, but in a rough upper limit estimation, the momentum 
for making this dip velocity gradient, P gra d, should be less than the product of the velocity 
difference and the whole mass of MMS1, P gra d < O.14M km s -1 . The momentum provided 
by the outflow is larger than P gra d; thus, this outer gas motion could be excited by the 
interactions with the outflow. 

The other point of discussion is the time scale before induced star formation is trig- 
gered in MMS1. One example of outflow-triggered star formation is L1551 NE in Taurus 



(lYokogawa et al.ll2003l ). In the case of MMS1, before star formation is triggered, it is neces- 
sary that a rise in pressure outside of MMS1 by colliding with the outflow reaches the center 
of the core. The velocity of this propagation equals the sound speed, or the shock speed if 
a shock has been excited. The sound speed under 10 K is 0.19 km s _1 . With this velocity, 
the signal of the collision arrives at the center in 7.2 x 10 4 yr. This is about an order of 
magnitude longer than the dynamical timescale of the outflow. Under this condition, trig- 
gered star formation would not have occurred yet. Next, we also consider the case in which 
a shock exists. In order to derive the shock speed, we need to know the density structure 
around MMS1, because the shoc k speed depend s on th e difference between inner and outer 



density of MMS1. According to iReipurth et al.l (119961 ). the difference between the 1.3 mm 



continuum intensity at the center and outside of the MMS1 is about a factor of 2. Under 
optically thin conditiona, the difference of the intensity is identical to that of density, and 
the shock speed equals half of the outflow speed, 5.3 km s _1 . With this speed, the shock 
would arrive at the center in 2.8 x 10 3 yr from the point of the collision. Considering the 
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time to reach the boundary of MMS1 from its departure from IRS4, the time for the outflow 
to arrive at the center of MMS1 is 8.5 x 10 3 yr, about the same as the dynamical timescale 
of the outflow. If the estimate of density structure is correct, compression has just reached 
to the center of MMS1 and collapse would just be triggered. Considering that we could not 
detect any infall and outflow signature around the source, it is possible that MMS1 is in a 
very young stage of its evolution. 

5. SUMMARY 

We observed an YSO condensation in the Cederblad 110 region of the Chamaeleon I 
dark cloud by CO(J=3-2), HCO+(J=4-3), H 13 CO+( J=4-3), CH 3 OH( J=7 K -Q K ), SO(J=8 9 - 
7s), and SiO(J=8-7) submillimeter lines with the ASTE 10m telescope. Results are as 
follows. 

1) We found a dense HCO + clump. This clump is an aggregate made up of three 
components associated with IRS4, IRS11, and Cha-MMSl. Although these sub-components 
do not appear in the integrated intensity map, we can distinguish them in the channel map 
and the position-velocity diagrams. 

2) We detected CO (J = 3 — 2) outflow components around IRS4 and IRS6, but not 
around Cha-MMSl, IRS11, or NIR89. Our observations with high angular resolution show 
the outflows previously known in this region are launched by IRS4, not by MMS1. This 
result suggests MMS1 is not a protostar but still in the prestellar phase. 

3) Each HCO+(J=4-3) line profile toward IRS4, IRS11, and MMS1 all have "red" 
asymmetry instead of blue asymmetry, which is a marker of infalling gas. However, our result 
does not necessarily mean there is no infall in these young sources, because the outflows from 
IRS4 and IRS6 are direct evidences of gas infall. 

4) The blue-shifted outflow from IRS4 seems to collide against MMS1. The outflow 
veers and the velocity declines just in front of MMS1. Though emission lines from shock 
tracer molecules such as CH 3 OH, SO, and SiO are not detected, it is plausible because the 
outflow from IRS4 is rather weak for excitation of these molecules. 

5) Momentum calculations show the outflow from IRS4 is strong enough to affect gas 
motion in MMS1, although the outflow is so young that the star formation process has not 
been triggered yet, or has just been triggered. This is consistent with the very young status 
of MMS1 and the lack of any outflows, jets, and masers. 
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Table 1. YSOs in Cederblad 110 Region 



Name 


R.A. 


(J2000) 


Dec. (J2000) IR class 


2~bol 


Lbdl 


References 


Cha-MMSl 


11 


06 


31.7 


-77 


23 


32 





20 


0.46 


1 



















36 


0.38 


2 


CedllO IRS4 


11 


06 


47.1 


-77 


22 


34 


I 


72 


1.0 


1 
















0/1 


59 


1.0 


2 


CedllO IRS6a 


11 


07 


9.8 


-77 


23 


05 


I 


260 


0.8 


1 


CedllO IRS11 


11 


06 


58.8 


-77 


22 


50 


I 




0.17 


3 


NIR89 


11 


06 


53.8 


-77 


24 


00 


I 








ISO-Cha I 97 


11 


07 


17.0 


-77 


23 


08 


II 








ISO-Cha I 101 


11 


07 


21.7 


-77 


22 


12 


II 








CedllO IRS2 


11 


06 


16.8 


-77 


21 


55 


III 


2700 


3.4 


1 



References. - (1) iLehtinen etaP toQlh : (2) iFroebrichl J2005h : (3) IPersi et al 
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Table 2. CO(J=3-2) Observations on YSOs 





Tub 


/ T MB dv 




(K) 


(K km s" 1 ) 


Cha MMS1 


6.94 ± 0.21 


14.4 ± 0.1 


CedllO IRS4 


8.39 ± 0.43 


34.7 ± 0.1 


CedllO IRS6 


11.37 ± 0.20 


24.3 ± 0.1 


CedllO IRS11 


6.50 ± 0.47 


15.5 ± 0.1 


NIR89 


7.86 ± 0.41 


17.6 ± 0.1 


ISO-Chal 97 


8.22 ± 0.21 


16.0 ± 0.1 


ISO-Chal 101 


6.94 ± 0.27 


14.1 ± 0.1 


CedllO IRS2 


7.80 ± 0.38 


14.9 ± 0.1 
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Table 3. HCO+(J=4-3) Observations on YSOs 





J- MB 


^LSR 


L\ V 


J J- MB Of 




(K) 


(km s _1 ) 


(km s _1 ) 


(K km s- 1 ) 


Cha MMS1 


1.64 ± 0.09 


4.49 ± 0.12 


0.68 ± 0.14 


1.11 ± 0.01 


CedllO IRS4 


3.00 ± 0.09 


4.46 ± 0.13 


1.55 ± 0.15 


4.58 ± 0.01 


CedllO IRS6 


0.65 ± 0.09 


4.62 ± 0.18 


1.53 ± 0.28 


0.73 ± 0.01 


CedllO IRS 11 


1.47 ± 0.08 


4.67 ± 0.13 


0.90 ± 0.15 


1.35 ± 0.01 


NIR89 


< 0.16 






< 0.35 


ISO-Chal 97 


< 0.17 






< 0.70 


ISO-Chal 101 


< 0.36 






< 0.54 


CedllO IRS2 


< 0.29 






< 0.30 



a ^LSR at the peak of the best-fit Gaussian profile. 
b FWHM of the best-fit Gaussian profile. 



Table 4. H 13 CO+( J=4-3) Observations on YSOs 





Tub 


^LSR a 




/ T MB dv 




(K) 


(km s _1 ) 


(km s _1 ) 


(K km s- 1 ) 


Cha MMS1 


0.26 ± 0.03 


4.63 ± 0.13 


0.36 ± 0.16 


0.11 ± 0.01 


CedllO IRS4 


0.13 ± 0.03 


4.53 ± 0.18 


1.72 ± 0.29 


0.18 ± 0.01 



a ^LSR at the peak of the best-fit Gaussian profile. 
b FWHM of the best-fit Gaussian profile. 
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Table 5. H2 Column Density and HCO + Abundance toward IRS4 and MMS1 





Cbeam 
^1.3mm 

(mJy) 


T D iV(H 2 ) 
(K) (cm" 2 ) 


X(HCO+) 


Cha MMS1 
CedllO IRS4 


370 (1 ) 
101( 2 ) 


20® 2.2 x 10 23 
21® 2.8 x 10 22 


2.6 x 10~ 10 
3.4 x 10~ 9 



References. — (1) 
Henning et all Jl993h : (3) 



Reipurth et al 



Lehtinen et al 



JIM; (2) 

J200lh 



Table 6. Parameters for Outflows from Ced 110 IRS4 and IRS6 





-Koutflow 


M 




li/f 




p 

1 








P 

1 






(PC) 


(Mr,) 


( vear") 




) 


fMokm s 


) 


(Mr-, fkm s" 


-1\2\ 
} ) 


(M^ km s _1 


y j 


IRS4 
























Blue 


6.5xl(T 2 


l.lxl0~ 2 


l.OxlO 4 


1.1x10 


-6 


7.1x10- 


-2 


2.4x10" 


-1 


7.1 xlO" 


-6 


Red 


2.5xl(T 2 


9.8xl0- 4 


2.9xl0 3 


3.4x10- 


-7 


7.7x10' 


-3 


3.1x10" 


-2 


2.7 xlO" 


-6 


IRS6 
























Blue 


1.6xl0- 2 


3.3 x 10~ 4 


7.8xl0 3 


4.3x10- 


-8 


6.5x10' 


-4 


6.7x10" 


-4 


8.4 xlO" 


-8 


Red 


3.1xl0~ 2 


5.7 x 10~ 4 


1.3xl0 4 


4.3x10- 


-8 


1.3x10" 


-3 


1.5x10" 


-3 


9.9 xlO" 


-8 
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Fig. 1. — Distribution of high velocity component of CO(J=3-2) line. The blueshifted (- 
4.0 km s- 1 < V LS r < 3.0 km s" 1 ) and the redshifted (6.3 km s" 1 < V hSR < 10.3 km s" 1 ) 
components are drawn with solid and dotted lines, respectively. The first contour level and 
contour interval are 1.53 K km s _1 for the blueshifted components and 1. 1 K km s" 1 for 
redsh ifted components. Cha-MMSl observed with a 1.3 mm continuum ( IReipurth et al. 



19961 ) is also shown with thick lines. The first contour level is 200 mJy, and the interval is 25 
mJy Small crosses indicate observed points. Each YSO is shown by filled circles (pre-stellar 
core and class I) or triangles (class II or III). The diamonds with the central dot denote the 
observed points of SO and SiO lines. The HPBW is indicated with an open circle. 
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Fig. 2. — Velocity channel maps of CO(J=3-2) in the outflow range. The central LSR 
velocity for each integral in units of km s -1 is shown at the top left corner of each panel. 
The first contour level and contour interval are 0.8 K. The markers are the same as in FigJTJ 
The gray region represents local minimum. 
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Fig. 3. — Integrated intensity map of CedllO in HCO + ( J=4-3) and the distribution of YSOs. 
The integrated velocity range is 2.0 km s _1 < Vlsr < 7.0 km s _1 . The first contour level 
and contour interval are 0.8 K km s" 1 . Dashed line shows the " absorption- free" integrated 
intensity contour (See §3.3.2). The open circle shows the HPBW. Filled circles and triangles 
indicate the positions of YSOs as in Fig{TJ The straight lines are also shown along which 
the position- velocity diagrams in FigJS] and [BJ are taken. 
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Fig. 4. — Velocity channel maps of HCO + (J=4-3) transition. The central LSR velocity for 
each integral in units of km s _1 is shown at the top right corner of each panel. The first 
contour level and contour interval are 0.5 K. The positions of YSOs are indicated with the 
same symbols as in FigHJ The gray region represents local minimum. 
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Fig. 5. — Position velocity diagram along the line (a) in Fig. [31 CO( J— 3-2) and HCO + ( J=4- 
3) are shown with thin and thick line, respectively. The first contour level and contour 
interval are 1.5 K for CO and 0.9 K for HCO + . The positions of IRS4, IRS6, and IRS11 are 
indicated with dashed lines. 
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Fig. 6. — Same as Fig. [U but for line (b) in Fig. [3j 
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Fig. 7.— HCO+(J=4-3) (multiplied by a factor of 2), CO(J=3-2) (with 2 K offset), and 
H 13 CO+(J=4-3) (only in IRS4 and MMS1 panels with 12 K offset and multiplied by a factor 
of 15) line profiles toward Cha-MMSl and YSOs in CedllO region. In panel a, the spectrum 
toward the edge of the map is indicated with a dotted line, and the borders between outflow 
and line core are shown with two vertical lines. We identified Vlsr < 3.0 km s" 1 and Vlsr > 
6.3 km s^ 1 as outflow. 
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